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Abstract 

We have computed persistent currents in a disordered mesoscopic ring in the 
presence of small electron-electron interactions, treated in first order pertur- 
bation theory. We have investigated both a contact (Hubbard) and a nearest 
neighbour interaction in ID and 3D. Our results show that a repulsive Hub- 
bard interaction produces a paramagnetic contribution to the average current 
(whatever the dimension) and increases the value of the typical current. On 
the other hand, a nearest neighbour repulsive interaction results in a diamag- 
netic contribution in ID and paramagnetic one in 3D, and tends to decrease 
the value of the typical current in any dimension. Our study is based on 
numerical simulations on the Anderson model and is justified analytically in 
the presence of very weak disorder. We have also investigated the influence of 
the amount of disorder and of the statistical (canonical or grand-canonical) 
ensemble. 
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Since the experimental discovery |T|-§| of persistent currents in phase-coherent mesoscopic 
rings pierced by a magnetic flux, a great deal of theoretical work has been devoted to the 
understanding of the order of magnitude of these currents, which is difficult to explain within 
a single electron theory. Both the amplitude of the average current measured in a many- 
rings experiment [Q, and the typical current measured in single ring experiments are 
found to be, at least in the diffusive regime, between 1 and 2 orders of magnitude larger 
than their expectations In order to explain this dicrepancy between experiment and 

theory, several authors have studied the influence of electron-electron (e-e) interactions on 



the persistent currents. Since the work of Altshuler and Aronov [|10| based on diagrammatic 
theory, it has been indeed known that transport properties of quantum coherent disordered 
systems are affected by e-e interactions. Inspired by these results, Ambegaokar and Eckern 
||TT| have shown that the presence of short-range interactions gives rise to a contribution 
to the average persistent current, whose sign depends on the sign of the interactions. The 
low flux current is found to be paramagnetic for repulsive and diamagnetic for attractive 
interactions. It has also been suggested that the typical value of the current could be 
strongly enhanced by interactions and disorder independent [^,|T^. In contrast with these 
results which are valid for 3D disordered rings in the diffusive regime, numerical ID exact 
diagonalisations have found that the current is decreased in the presence of repulsive nearest 



neighbour interactions The case of long range Coulomb interactions has also been 



investigated ||T5|,|T6[ , with the result that the current is decreased for a very weakly disordered 
ring and increased at higher disorder. However from these studies performed in ID it is very 
difficult to anticipate the behavior of a realistic 3D multichannel ring, which is very difficult 
to handle exactly numerically. 

These considerations have motivated the present work where we have studied the influ- 
ence of a small e-e interactions, on disordered rings in ID and 3D. We have compared the 
cases of a pure contact Hubbard and a flrst neighbour interaction. The e-e interaction is 
treated within flrst order of perturbation on the Anderson model as follows: 
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^ = H [tijCfCj + C.C.] + WiCfCi (1) 

with tij = exp (2?7r$/$o(a;j — Xj)/Lx), where i and j are nearest neighbours, and where 
Wi is the onsite disorder potential uniformly distributed in the interval [—W/2, W/2]. This 
hamiltonian is exactly diagonalized in the presence of disorder . 

We first consider the one-band Hubbard model, characteristic of contact interactions 
described by the hamiltonian Hint = ^ X! ^jT'^ji with njo- = c^Qo- and f/ > for repulsive 

i 

interactions. The first order correction to the energy of a N electrons state reads : 

p 

5E{2p) = ^ ii N = 2p 

p (2) 
SE{2p + 1) = 6E{2p) + J2 W(«> P + 1) if iV = 2p + 1 

a=l 

Here U{a, (3) describes the interaction energy between electrons in the states a and (3 which 
are single electron eigenstates of the Anderson hamiltonian. 

W(a,/5) = f/El^«WlV/3»l' (3) 

i 

where ipai^) denotes the value of the one-electron wavefunction at the site i. We can also 
write the correction to the ground state with 2p electrons in the form: 

5E{2p) = UY.nl{^ (4) 



where np{i) = ^ is the electronic density at the site i with the first p eigenstates 

occupied by only one electron. Expression (|^) is very similar to the result obtained by 
Argaman and Imry We deduce the correction to the current given by: 

and compute the canonical (C) and grand-canonical (GC) average of this current in 3D in the 
diffusive regime (the averaging is done over the filling, between 1/4 and 3/4 of the spectrum). 
In the GC case, only configurations with an even number of electrons contribute. In fig. ^a 
these quantities are compared to the average C and GC currents without interactions. We 



can see that, unlike the non-interacting case, both statistics give identical results for (SI) 
which has the same paramagnetic sign (in zero flux) as the canonically average current in the 
absence of interactions. Furthermore, the effect of interactions is to enhance the current, 
which is in agreement with analytical results From the first order correction to the 

energy , it is also possible to compute the quantity {6{P{N))) = 2{I{N)6I{N)) which is the 
first order correction in U to the typical current (/^). We observe in fig. ||b that (ISI), in 



both C and GC statistics, is always positive and of the order of y{P){6P), which means 
that there is a strong correlation between the current of N electrons and the correction to 
this current for the same number of electrons. We have made this calculation for different 
values of the disorder. The first harmonics of average and typical currents are presented in 



fig. D as a function of the disorder W. We observe that (J2), {SI2), y {If) and y 
decrease as oc le, the mean free path. This dependence of {6I2) with the disorder is in 
complete agreement with analytical calculations Nevertheless, as already pointed out 



IP, the disorder dependence of the C average {I2) for this small ring disagrees with analytical 
predictions [§-0]. In conclusion we have found that the effect of a contact interaction results 
in an increase of both average and typical currents. In the following we demonstrate this 
result analytically in the weak disorder limit. Eigenstates of the Anderson hamiltonian in 
ID without disorder are plane waves of the type ipn,<!>ix) = exp(i2Tc(n — ^/^q)x/L)/^/L. In 
the presence of a very weak disorder, the degeneracy occuring at $ = and $ = $o/2 is 
lifted, resulting in the following wavefunctions: 



^n,M = j2/Lcos[27r{n-^/^o)x/L + Xn] 

'^n,M = ^j2/Lsm[27i{n-<^/%)x/L + Xn] 
where Xn is a phase factor which depends on the particular realisation of the disorder. We 

can easily deduce the following value of 5E{2p) : 

(p^ + 1) /L for ($ = and p = 2N) or ($ = $o/2 and p = 2N + 1) 
SE{2p) = J V F ; V u/ F ^ (7) 

p'^/L otherwise 

The resulting presents an antisymmetric paramagnetic peak in the vicinity of $ = 

(resp. $0/2) for an even (resp. odd) number of electron pairs, in agreement with the 



symmetry property of the ID electron spectrum, according to which: I{N, $) ~ I{N + 
1, $ + $0/2), vahd in the hmit ^ 1. This behavior gives rise to a current contribution SI 
which has the same sign as in the absence of interactions (i.e. paramagnetic in zero flux for 
2p electrons) in agreement with our numerical findings (see fig. ^). 

In comparison with the Hubbard model we also considered the short-range nearest neigh- 



bours interactions studied by Bouzerar et al. The corresponding hamiltonian reads : 



'Hint = cfcicf_^iCi^i. Two terms appear now in the first order correction to the energy : 

i 

a direct and an exchange term, so that the final correction reads : 

SEi2p) = u(2j2npii)np{z + l)-J2 E C Wa(^ + 1)^;(^ + 1) ) (8) 

\ i i a,/3=l / 



Note that unlike the case of ref. [T^, here spin has been taken into account. However we 



have found that the exchange term is always much smaller than the direct one, which means 
that the expression obtained for spinless electrons would only differ by a factor two from 
what we have calculated. 

In ID, we find a current contribution which, for each value of the number of electron 
pairs (which varies from to |L, L being the length of the ring) is opposite in sign to 
the current in the absence of interaction, in agreement with ref. [0]. The resulting average 
current gives rise to a diamagnetic contribution (see fig. |b). 

On the other hand in 3D, the same type of interaction gives rise to an average paramag- 
netic contribution to the current, which is shown in fig. ^ for different values of the disorder. 
We observe like previously that the current does not depend on the statistics considered, 
but its amplitude is much larger than the value obtained for a contact interaction. The 
dependence of the average current on disorder is also very weak in contrast with the results 
obtained for the Hubbard interaction (see fig. 0a). If we now study the results obtained 
for the corrections to the typical current in 3D we see that for most values of the fiux and 
for ballistic and diffusive regimes (ISI) is negative for the C statistics except in the vicinity 
of $ = and $ = $o/2 where it is slightly positive. In the GC ensemble (ISI) is always 
negative (see for instance fig. § obtained in the diffusive regime) . The typical current is thus 



diminished in the presence of first neighbour interactions both in ID and 3D. The variation 
of the first harmonics versus disorder are presented in fig. |^. We observe that (5/2) generally 
does not vary very much with disorder. 

We have also performed an analytical computation of 5E{2p) (as previously discussed 
in the Hubbard case) for a ID ring in the limit of very weak disorder, and obtained for the 
first order energy correction : 

[flip) + f2{p)]/L for ($ = and p = 2N) or ($ = $0/2 and p = 2N + 1) 
6E{2p) = { (9) 

fi{p)/L otherwise 



with 



flip) = 2p'-( 'J^PI^y and f,{p) = \ cos2pn/L (10) 

Note the strong dependence of SE{p) on the filling : the sign of cos2p7r/L determines the 
flux dependence of SE{2p) and consequently the sign of SI{N, $). When p varies between 

and |L, cos2p7r/L is negative. It results from this that 5/($) presents an antisymmetric 
diamagnetic peak near $ = ($o/2) for an even (odd) number of electron pairs. In the 
weak disorder limit this calculation can be generalised to a multichannel ring. In that 
case we can separate the contribution of the interactions to the persistent current into its 
longitudinal and transverse part. The longitudinal component can, depending on the filling, 
lead to diamagnetism; however, the transverse part of the interaction always gives rise to 
paramagnetism, since it is identical to the Hubbard contribution. This explains why the total 
contribution is paramagnetic. Nevertheless, note that the case of anisotropic interactions 
may lead to the opposite conclusion. So we have obtained a confirmation of our numerical 
results which show that a first neighbour interaction gives a diamagnetic contribution to the 
average current in ID, and also tends to decrease the value of the typical current. Only this 
last result remains valid in 3D, where the average current is paramagnetic. 

From our perturbative study of the influence of e-e interactions on the persitent current in 
a disordered ring, we have shown that a repulsive Hubbard interaction enhances the typical 
value of the persistent current in any dimension. It also contributes to the average current 
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with a paramagnetic sign in zero field. These extra both typical and average contributions 
to the current decrease with disorder just like the single electron quantities. There is thus no 
evidence that the Hubbard interaction could attenuate the sensitivity of the currents upon 
disorder. 

A nearest neighbour interaction, on the other hand, tends to decrease the value of the 
typical current, while the sign of the average contribution depends however on the dimen- 
sionality of the ring : it is diamagnetic in ID and paramagnetic in a multichannel ring, with 
very small disorder dependence. These differences between Hubbard and nearest neighbour 
interactions in the one dimensional limit have also been found analytically. 

This work has strongly benefited from the help and suggestions of A. Altland, Y. Gefen, 
T. Giamarchi and G. Montambaux. We thank also J. A. Cowen for having revisited the 
manuscript. Numerical simulations have been performed using CRAY facilities at IDRIS 
(Orsay). This work was partly supported by grant from DRET No. 92/181. 
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FIGURES 

FIG. 1. (a): flux dependence of C and GC corrections to the average current, compared to 
their values without interactions, (b): C and GC corrections to the typical current compared to 
the typical C current without interactions. All curves are calculated for a 30 x 4 x 4 ring, for a 
contact (Hubbard) interaction with W = 2 and U = 1. 

FIG. 2. Dependence of different harmonics of the correction to average or typical current on 
disorder. The effect of Hubbard (subscript "H") or nearest neighbour (subscript "NN") interactions 
are compared with the average and typical currents without interaction (all of them for a 30 x 4 x 4 
ring with U = 1). (a): second harmonics of the average current with Hubbard, nearest neighbour 
or no interaction, (b): first harmonics of the typical current with or without Hubbard interactions, 
(c): first harmonics of the typical current with or without nearest neighbour interactions. 

FIG. 3. Flux dependence of 61 for a Hubbard (a) or nearest neighbour (b) interactions for 
all the even fillings of a ID ring of length L = 16, from to |L. The circles correspond to the 
average within the same range 

FIG. 4. (a): flux dependence of (SI) for a nearest neighbour interaction for three amplitudes 
of disorder, in the C ensemble in 3D. (b): flux dependence of (negative) C and GC corrections to 
the typical current compared to the typical C value, for a 30 x 4 x 4 ring with nearest neighbour 
interactions, with W = 3 and U = 1. 
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